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We study hard exclusive electroproduction of two pions in the QCD factorization approach at 
next-to-leading order (NLO) in the strong coupling. The pion pair can be produced both in an 
isovector and in an isoscalar state. The angular distribution of the produced pion pairs allows 
one to project out a component that depends on the interference of the isovector and isoscalar 
channels. Using specific models for the involved generalized parton distributions and two-pion 
distribution amplitudes we investigate the angular distributions of the pion pair in NLO and com- 
pare them with HERMES data. The differences between the LO and NLO results are moderate 
and the agreement with data is satisfactory, though not perfect. Expecting new results from COM- 
PASS collaboration we perform the calculation also for the COMPASS kinematics. 
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1. Introduction 

In the last years it has been realized that a large class of hard exclusive reactions can be treated 
in QCD factorization framework, absorbing all non-perturbative soft physics in suitable generalized 
parton distributions (GPDs) and distribution amplitudes (DAs). GPDs encode valuable information 
about hadron structure, including some aspects which cannot be deduced directly from experiment, 
like the transverse spatial distribution of partons and the total angular momentum. For more details 
we refer to the reviews [jl], ^, 

The large amount of information contained in GPDs implies that much and diverse data is 
needed to determine their functional form in the three variables on which they depend. One of the 
channels in which data is already available is the exclusive electroproduction of pion pairs This 
was already studied by some of us ^ in the past in LO and is now analyzed in NLO. 



2. The amplitude of the process 

The analysis of the di-pion electroproduction is reduced essentially to that one for the subpro- 

cess 

t{q) +N{p) ^ n'^in + %\k^) +N{p'), (2.1) 

where a longitudinally polarized virtual photon with momentum q hits a nucleon N with mo- 
mentum p and produces a final state nucleon N with momentum p' and two pions with momenta 
k" and k^. We use the conventional variables 

q^ = -Q\ A = p'-p, A^=t, ^Bj = x^, (2.2) 

2p-q 

and denote by mj^,mN and mj^n the pion, the nucleon and the di-pion mass, = {k" + ^^)^. We 
consider the limit that virtuality of the photon is large, 

> |f|,m^,m^,m^^ , (2.3) 

in fact, is not larger than the energy scale W^. In this case the amplitude, according to the 
factorization theorem in may be written as a convolution of hard coefficient function and soft 
parts, parameterized by GPDs, and 27rDAs. Note that the subprocess with longitudinally polarized 
photon gives the leading contribution to the reaction, the contribution of transverse polarization is 
suppressed by a power \/Q. 

At leading twist the pion pair can be produced both in isoscalar and isovector states 

jn+n- ^j,I=0 ^jl=\ ^ (-2.4) 



where 

1 

eW(z,T)f^W(T,§,0 (2.5) 
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Figure 1: Typical NLO diagrams for di-pion production in the isovector (upper line) and isoscalar (lower 
line) state. The diagrams related to each other by crossing are displayed one upon another. 



^,^0 ^ W4^ ]dr]dzl^ efFf(-Hr,^,t) [G(-)(z,T)<I>f «(z, 



NcQ^ 



-1 



f=u,d 



f'=u,d,s 



ntj, 



Here t" = 1 and t| = — 1 for the up- and down-quark, ^ = xbj /(2 — XBj ) is the skewness parameter, 
Nc the number of colors and a the fine structure constant. Q^^\ G^^^ and R^"^^ denote hard coef- 



ficient functions, whereas F and <t> stand for GPDs and 27rDAs. The integration in (2.5) runs over 
parton momentum fractions z and z- Due to charge conjugation invariance the production of an 
isovector pion pair, described by a C odd quark DA, is mediated by gluon and C even quark GPDs, 
and F^^\ In contrast, the C odd quark GPD F'^) contributes for the production of an isoscalar 
pion pair, which is produced by hadronization of a a gluon pair or a C even quark pair. 

The coefficient functions represent the amplitudes for the scattering of coUinear partons. Typ- 
ical NLO diagrams are shown in Fig. [I], where the graphs in the first, the second and the third 
columns contribute to R^^\g^^^ and Q^^^ respectively. 2'-+' were calculated in the MS 

scheme in [pi], where electroproduction of light vector mesons was studied at NLO. For example 



2y+l 

y(>'+i) 



Ri+){z,r) = aj{^l)CF^£g[z,'-^) , (2.6) 
Zin2(_3;) _ ^ln2(3;+ 1) + (in (§) - l) iyln{-y) - (3;+ l)ln(y + 1))" 



V{z,y) ^ .yln(-y)+(.y+l)ln(y+l) , y(.v+l)+(>-+z) 



y+z 



-Wiz,y)} + {z^-z}, 
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Figure 2: Definition of the polar angle d in the cm. of the pion pair. 

where }Xr {}Xf) is the normalization (factorization) scale, Cp = {N^ — \)/2Nc, and 

V{z,y) = z\n{-y) +zln(3; + 1) +zln(z) +zln(z) 
W{z,y) = Li2(j+ 1) -Li2(-y) +Li2(z) -Li2(z) +ln(-3;)ln(z) -ln(3; + l)ln(z) . 

The hard amplitudes for pion production in the isoscalar state can be obtained by crossing from 
those for the isovector state. They correspond to each other after the interchange of the f-channel 
and the ^-channel partonic pairs and the appropriate interchange of the relative partonic momentum 
fractions. One can convince oneself of this relationship easily by comparing the typical NLO 
diagrams on the upper and lower hne in the Fig. |l|. The prescription for the interchange of the 
momentum fractions reads 

z-^, z-^, (2.7) 
then e(-) (z, t) = 2^+) (^^, (2z - 1)) , and similar relations hold for R'^-\G^-\ 

3. Numerical results 

We will present results for normalized Legendre moments (NLMs) defined as a convolution of 
differential cross section 

/dcos0 P„{coid)'^''"""'' 



dcos B 

{Pn{cOSd)f'''" = (3.1) 

/dcose ^ 

J dcos 9 

-1 

with a Legendre polynomial P„(cos 6), where the polar angle is defined in Fig. |2[ 

The GPDs in our calculation are modeled using the ansatz for double distributions suggested 
by Radyushkin |^. Our models for 27rDAs reads 

^=\z,^,m^^) = 6z(l -z)(2C - l)F^im^^) , (3.2) 
0^=0 ^ , ) = _ H2^2( i-,)(2z-l)[^foMPo (cos d)-^f2{m^M<^osd)\, 

^^{z,^,m^^) =-60M^zHl-zf ?-^Mm^Mcosd)-^f2{mnMcose)\ , 



where 2^ — 1 = j8cos6, the pion velocity in the di-pion cm. is p = y^l —4m^/m^^, and Uf is 
the number of active flavors. and M2 represent the momentum fractions carried by quarks and 
gluons in the pion target. The isovector pair is produced in a P-wave, and F-,j^{mjt„) is the timelike 
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electromagnetic pion form factor. The isoscalar pair may be in an 5- or D- waves, and foim,^^) and 
fii^nn) are the corresponding Omnes functions for S- and D- waves. Omnes functions develop an 
imaginary part above two-pion threshold. In the region where pion-pion scattering is elastic, the 
phases of fQ{mnn) and /2(?M;r;r) coincide with the pion phase shifts 5q"°, 8{""^ (Watson theorem). 
For more details see [0, ll. For higher di -pion masses the phases of Omnes functions do not 
coincide with the pion phase shifts. In the present analysis we neglect inelasticity and use the pion 
phase shifts as an input in dispersion relations to reconstruct the Omnes functions. We use two 
sets (SI and S2) of parameterizations for the Omnes functions. In the set SI we calculate Omnes 
functions using dispersion relations with two subtractions and the fits [^] for the pion phase shifts. 



For the subtraction constant we used the result of [10|. The set S2 is the same as that one used in 

Odd Legendre moments are proportional to the product of isoscalar and isovector amplitudes 
(Pi3)oc yt/cos0A,3(cos0)/?e 1(7^=^7^=°} (3.3) 
and sensitive to the interference of P-wave with S- or D- wave, 

(A) - Re{F;{m^n){cih{mnn)+C2f2{mnn))] , {Pi) - Re{F;{m^n)f2{mnn)] ■ (3.4) 

These observables provide access to a small isoscalar amplitude. 

In Fig. ^ we compare the results of our calculation (with set S 1) with HERMES data for NLMs 
on hydrogen and deuterium targets [Q]. Note that the difference between the LO and NLO results 
is not big, which may indicate a fast convergence of the perturbative series. Our leading twist 
predictions for Pi are in reasonable agreement with the experimental data for the hydrogen and 
deuterium targets. Since such experiments in principle allow one to test the gluonic content of the 
nucleon, we plot in addition the corresponding results when we omit the two-gluon exchange in 
the f-channel. The effect of the gluon contribution is noticeable, hence improved data in the future 
may provide constraints on the gluonic content of the nucleon. 

In Fig. we present our results for kinematics typical of the COMPASS experiment, i.e. for 
larger and smaller xbj. We plot here the results obtained with different Omnes functions, with 
and without account of two-gluon exchange the f-channel. The predicted values of the NLMs are 
smaller for COMPASS kinematics than that ones for HERMES. Another observation is that the 
gluon GPD plays here even a more important role. 

To estimate the scale uncertainties of our NLO calculation we plot in Fig. ^ for HERMES 
kinematics the "longitudinal" combination of the moments, which projects out the state with van- 
ishing total helicity of the pion pair, calculated with set S 1 and different settings in the hard part 
for the factorization and the renormalization scales. We see that this scale uncertainty is not large. 



4. Summary 

We studied di-pion electroproduction in the QCD factorization approach at NLO. Numerical 
results obtained for HERMES and COMPASS kinematics show that NLO corrections (at least for 
normalized Legendre moments) are small. Hence with improved experimental data from such 
experiments it will be possible to put the constrains on both GPDs and 27rDAs. One can also get 
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Figure 3: The m,iji dependence of the normalized Legendre moments (Pi) (upper panels) and (P3) (lower 
panels) for hydrogen (left panels) and deuterium (right panels). The curves show our results at LO (dashed 
black), NLO (solid black) and at NLO without two-gluon exchange in the f-channel (dotted). The average 
kinematics for the HERMES data [|l is (xbj) = 0.16, {Q^) = 3.2 (3.3)GeV2, and (-r) = 0.43 (0.29)GeV2 
for hydrogen (deuterium). The histograms show the systematic uncertainty of the measurement. 
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Figure 4: The normalized Legendre moments calculated for a deuterium target in the kinematic region 
relevant for COMPASS ( xbj = 0.08, -f = 0.27 GeV^ and = TGeV^) with sets SI (solid) and S2 (dashed 
lines). The dotted and dashed-dotted lines represent the calculation without two-gluon exchange in the 
/-channel, for sets S 1 and S2 respectively. 
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Figure 5: The "longitudinal" NLM combination calculated for different renormalization and factorization 
scales, hr = Hf = H- 

additional insight into the gluonic content of the nucleon. To be on the safe side with the leading 
twist approach it is, however, important to go to larger Q^. 

In the future we want to improve our approach to Omnes functions by taking into account 
inelasticity effects. 
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